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X-ray Coherent diffraction interpreted through the fractional Fourier transform 
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Diffraction of coherent x-ray beams is treated through the Fractional Fourier transform. The 
transformation ahow us to deal with coherent diffraction experiments from the Fresnel to the Fraun- 
hofer regime. The analogy with the Huy gens- Fresnel theory is first discussed and a generalized 
uncertainty principle is introduced. 

PACS numbers: 



DIFFRACTION OF A RECTANGULAR 
APERTURE 

Almost two centuries ago, Fresnel has shown a beau- 
tiful consequence of the wave character of light. If an 
opaque disc is located in front of a monochromatic light 
source, a bright spot is observed at a given distance Zh at 
the center of the shadow of the disc[l]. This phenomenon 
is known as the Poisson bright spot. In the opposite sit- 
uation, the diffraction pattern of an aperture gives rise 
to a minimum of intensity if the detector is located at a 
very specific distance. This minimum of intensity will be 
called the dark spot in the following. 

The dark spot results from a destructive interference. 
It is located in the Fresnel regime where the beam re- 
mains almost parallel and where the amplitude oscillates 
rapidly. The diffraction pattern of a squared aperture 
simulated from the Fresnel's integral is displayed in Fig[T] 
where the dark spot is clearly visible. 

The Fraunhofer regime is reached for distances be- 
tween the detector and the diffracted object [2], larger 

than z » ^5f~* -"-^ ^^^^ regime, the diffracted amplitude 
is proportional to the Fourier Transform (FT) of the am- 
plitude at the diffracted object position: 
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with <7 = In the case of a rectangular function s(x) 
which is only non-zero in the interval [— f ,f ], we obtain: 



I{q) = F*{q)F{q)^[ 



(2) 



The diffracted intensity by a rectangular function is 
proportional to the well-known cardinal sinus function 
squared, in perfect agreement with measurements. The 
measurement of the diffraction pattern of a 2fimx2fim 
slit by a x-ray beam in the Fraunhofer regime is dis- 
played in Fig 2. In this regime, the amplitude varies 
very smoothly and the beam width A is proportional to 
2;: A oc \z/a. 

This old problem of optic has recently recovers a great 
interest by the scientific community using coherent X-ray 



diffraction. Indeed, coherent X-ray beams are obtained 
from weakly coherent synchrotron sources by cleaning 
and collimating the beam thanks to rectangular slits [3|. 
Concretely, two sets of rectangular slits are usually used 
along the optical path. The downstream one is located 
at few centimeters from the sample and is opened at few 
tens of micrometers. Within this setup, the knowledge of 
the dark spot position is important because it is usually 
located at few tens of centimeters downstream the slit. 
As a consequence, the sample may be located either in 
the Fresnel or in the Fraunhofer regime as respect to 
aperture and wavelength, which may strongly influence 
diffracted patterns. 
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FIG. 1: Diffraction of a rectangular aperture in the Fresnel 
regime and corresponding profiles simulated from the Fresnel 
integral (eq. [5| or from the Fractional Fouier Transform (eq. 
[4|). The dark spot is located at Zh and the minimum beam 
width at 

The dark spot position is not obvious to calculate ana- 
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FIG. 2: Diffraction pattern of a 2/j.m x 2/j.m squared aperture 
obtained at A = 1.5 A (Log scale). The camera is located 
2,22m downstream. A diagonally placed absorbing wire has 
been used across the direct beam to protect the camera (from 
reference Hi). 



lytically. In first approximation, Zm can be estimated by 
assuming that it corresponds to the distance z where the 
beam width in the Fraunhofer regime (A ^ A^) is equal 
to the beam width in the Fresnel regime (A ^ a), that 



A ' 



From numerical simulations, the dark spot position is 
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located close to z^ ^ and the location where the 
beam width is minimum at ^ fx (see Fig [1]) . 



FRACTIONAL FOURIER TRANSFORM 

We first discuss in this section the relation between the 
Huygens- Fresnel theory and the fractional Fourier trans- 
form, introduced by NamiasQ. The previous slit diffrac- 
tion in Fig [T] can be obtained from the fractional Fourier 
Transform. Within this framework, the resolution of the 
Fresnel's integral is nothing else but the resolution of the 
quantum harmonic oscillator. 



Fractional Fourier Transform: an operator 



If the same operator is applied two times, we obtain: 

^'im = /(-^)- 

The operator T has to be applied four times to recover 
the original function: 

as illustrated in Fig |3] from a 2D picture. Within this 
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FIG. 3: A classical Fourier transform is an angular transfor- 
mation of 7r/2. The operator has to be applied four times to 
recover the original portrait. The starting point of the simula- 
tion is a painting by Manet "Berthe Morisot" Musee d'Orsay, 
Paris (displayed in black and white and considering bright- 
ness variations as density variations). Only the intensity is 
represented in the reciprocal space. 

framework, the standard Fourier Transform (FT) is an 
angular transformation of 7r/2 in the (x,q) plane. A gen- 
eralized Fourier Transform can thus be developed for any 
angle of rotation a from zero to ^. The integral form of 
the fractional Fourier transform (in ID case) can be writ- 
ten asjil, 0: 



/oo 
K{x,u)f{u)dt 
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with 
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For a = 0, it can be shown that this expression is equiv- 
alent to the identity [7]. The standard Fourier transform 
(eq. d]) is easily obtained for a = 7r/2. 



The most natural way to introduce the Fractional 
Fourier Transform is to note that the Fourier Transform, 
defined in eq(TJ can be written as an operator T acting 
on function: 



The fractional Fourier transform and slit diffraction 

The Fractional Fourier Transform (FrFT) in 2D case, 
so with appropriate factor in front of integral is: 

-^a/lq) = 7, — ^exp(-^q^ cot a)... 
27r sm a 2 
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It has to be compared to the diffracted field amphtude. 
We consider a diffracting plane E (containing the slit) 
and we observe the field amplitude at distance z of E 
on a plane screen 11 orthogonal to the axis. The field 
amplitude at a position on 11 defined by the variable ^ is: 
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The similarity between equation |4] and [5] suggests to 
write diffraction in terms of FrFT. Papers by Pellat-Finet 
d, [qI give the relation between the two expressions. The 
method is based on the use of intermediate spherical sur- 
face on which the amplitude of the field is considered, 
then within conditions on curvature and position of these 
surfaces, field on it can be related by Fourier transform 
or fractional Fourier transform. 



The Fraunhofer regime 

The simplest example is the case of the Fraunhofer 
regime. Consider two spherical surfaces (see Fig |4^). 
The first one is the object surface S of radius z pass- 
ing through the origin Qs (the middle of the slit if the 
object is a slit) and with a center Qe on the optical axis 
at positive distance z of Qs- The other spherical surface 
a screen where the field is observed, contains Qe and 
has center ^l^, so his radius is — z if oriented from surface 
to center along the axis. Points on these surfaces are de- 
fined by the coordinates of their projection parallel to the 
axis on the planes orthogonal to the axis in and Qe- 
In this case the field Us on S and Ue on E are related 
by a Fourier transform: 



(6) 



Using reduced variables p = (|^)-^/^r and p' = (|^)^/^r' 
as arguments of scaled functions Vs{p) and Ve(pO- 



VE{p')=iT./2[Vs]{p') 



(7) 



The field on E is the Fourier transform of the field on S. 
If z ^ oo this gives the limit of the Fraunhofer regime: a 
plane object is Fourier transform into a plane diffraction 
at infinity. Spherical surfaces in place of planes have to 
be considered for finite z. The Fraunhofer regime corre- 
sponds to the large z domain in which phases shift be- 
tween planes and spherical surfaces remain small. 





FIG. 4: Fraunhofer diffraction (a): the field on E is the Fourier 
transform of the field on S. Fresnel and fractional Fourier 
transform (b): the field on S is change into the field on D 
by a fractional Fourier transform defined by a with cot (a) = 
(1 — /i)//i. The field on S can be transform into the field on E 
combining two fractional Fourier transform defined by angles 
a and /S such that a -\- /S — 7t/2. 



The Fresnel regime 

In the FrFT, a = corresponds to z = and a = 7r/2 
corresponds to 2; ^ oo. It is tempting to get a rela- 
tionship between z and a. Comparison of two terms 
exp(— ^^) in equation m and exp(— ^^^f'^ ) in equation [5] 

suggests 
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but other terms in these equations contradict this expres- 
sion mainly for small a, in the Fresnel regime. There is 
a difficulty comparing variable q which defines a wave 
vector in the reciprocal space with ^ in real space. 

The solution to this problem has been given by Pellat- 
Finet. In place of the spherical surface through Qs with 
radius z a surface C with a smaller radius pz also through 
is considered. Call its center. We consider also a 
spherical surface B through Qc and center l]^, so with 
oriented radius —pz. Field transfer from S to C implies a 
quadratic phase factor, and transfer from C to 5 is again 
a Fourier transform. Consider a reduced field function Vb 
of a reduced variable p^\ 



Vb{p") 



2tt cos(q;) + sin(a) 



(9) 



The reduced variable are defined by p = (|j)-^/^r and 
p" = (^)^/^r'^ Notice that variables r and r" give posi- 
tion of points on spherical surfaces C and B. Such points 
are defined by coordinates of their orthogonal projection 
on a plane orthogonal to the optical axis. 

The p parameter smaller than one, is related to a by 



cot(a) 
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cos(a) + sin(a) * 
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In order to have a fractional Fourier transform with the 
correct quadratic phase factor, an other spherical surface 
D is needed. A sphere through Clc with positive oriented 



4 



radius R such that = (/i^ + (1 - /i)^)^/(l - 2/i). Note 
that surface D corresponds to a plane when /i = ^ . Then 
we have: 

Vd{p") = e'«(cos(a) + sm{a))Tc[Vs]{p")- (11) 

It can be shown that transfer from the field on D to the 
field on using intermediated spherical surfaces D' and 
E' to have correct phase factors, is given by a fractional 
Fourier transform where j3 defined by the coefficient 
fi' = 1 — II IS j3 = / 2 — a. 

This prove the continuity of the transformation fol- 
lowing that of the fractional Fourier transform, = 
J^oc + ^(3 ' Notice that this continuity is false for the Fres- 
nel transformation. 



Relation with the quantum harmonic oscillator 

The most remarkable property of the fractional Fourier 
transform is that a possible choice for the eigenfunctions 
of the operator J^ct is given by the set of normalized 
Hermite-Gauss functions[5, To|, U], similar to the orthog- 
onal harmonic oscillator basis: 
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where Hn{x) = (-1)"^ exp[x^/2]^ exp[-x 
eigenvalues of J^cx are e*^": 



72]. 
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Extension to two dimensions is straight forward. Obvi- 
ously, J^aJ^p^n{x) = Ta+(3^n{x). 

Consequently the diffraction problem can be mapped 
on that of the harmonic oscillator. The Fresnel's inte- 
gral or the fractional Fourier transform applied to the 
diffraction of a squared aperture is equivalent to the res- 
olution of the Schrodinger equation with time in a har- 
monic potential, with a rectangular function as initial 
time condition [12|. Then time evolution is related to the 
a variation. In terms of quantum mechanics, the frac- 
tional Fourier transform offers the way to continuously 
switch from the position space to the impulsion space. 



Gaussian beams 

It is interesting to consider the behavior of a Gaus- 
sian beam propagating from the Fresnel to the Fraun- 
hofer regime (see Fig [5^). Since the first eigenfunction 
of the basis is a Gaussian, it could be expected that a 
Gaussian beam will remain Gaussian if its width corre- 
spond to that of this eigenfunction. Following the Pellat- 
Finet approach, we have introduced reduced variables 
p = (||)^/^r. It is with such variables that the width of 
a beam have to be given. Then we have to consider, as 



given in Fig HI a Gaussian beam modulating a spherical 
wave S with center ^l^; at distance z. If the Gaussian 
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amplitude on S is proportional to exp(— with a = 1 
like for the first eigenvalue of the harmonic oscillator, this 
width will be kept propagating from S to E through D. 
Expressed using r, or r^' variables the width is (|^)^/^. 
With numerical values A = lA and z = 1 m, the constant 
width is 4 /im. The physical meaning of this value is that 
a beam (A = lA) falling on a curved mirror designed to 
focused at Im, with a width of 1/im just after the mirror, 
have in fact a constant width. 

What is the behavior of narrower or wider beams 
(cr 7^1) than the first eigenfunction ? The 1-D fractional 
Fourier transform f{x) of a Gaussian function /(x). 
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is: 



:exp(- 
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7ri/4y^(sina - ia'^ cos a) 2 a'^cota + i 

(13) 

Note that the amplitude Fa{q) remains a real function 
and is invariant with a when <j = 1. In the general case, 
it appears that the transform of a narrow function is wide 
and reverse, but it appears also that the Fa{q) function is 
not a Gaussian for a G]0,7r/2[. Nevertheless the module 
of this function remains a Gaussian in agreement with 
the experiment. The amplitude of the beam has a width 
6= ^ Vl - cos(2a) + cr4(l + cos(2a)). 



Fractional Fourier transform and the uncertainty 
principle 

In terms of quantum mechanics, every func- 
tion f{x) with a normalized probability density 
function (J^^ \f{x)\'^dx = 1) and its Fourier transform 
F{q) fulfilled the inequality: 



Var[/(x)] X Var[F((7)] > 



1 

4' 



(14) 



with Var[/(x)] = /^(x - x2)2|/(x)|2dx and 
X = J^^x\f{x)\'^dx. This principle is a direct property 
of the standard Fourier transform. As respect to the 
FrFT, the uncertainty principle appears to be a peculiar 
case for a = 1^ and can be generalized for any a: 

Var[/(x)] X Var[J'"/(x)] a e [0,7r/2] (15) 



The generalized uncertainty principle for a Gaussian beam 

Let's first consider the case of a Gaussian beam as 
given in equation [T2l which fulfilled J^^\f{x)\'^dx = 1 
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FIG. 5: a) Propagation of a Gaussian beam treated by the 
FrFT for a G [0, f ]. b) Product of variances Var[J^a/(x)] x 
Var[J^/3/(x)] with {a, 13) G [0, f ]. c) Coherent diffraction of 
a rectangular function treated by the FrFT as respect to 
a G [0, 1^]. The minimum of the rms variance is indicated 
by an arrow, d) Product Var[J^c^/(x)] x Var[J^/3/(x)] with 
G [0, f ] if the variance is measured over a hmited area 
centered at the beam maximum. The minimum product value 
is indicated by an arrow. 

2 

and Var[/(x)] = The Fourier transform F(q) gives: 

F{q)=:F.f{x) = ^eM-^) (16) 
7r4 2 

with Var[F(g)] = The minimum value is obtained 

for the Gaussian probability function: 

Var[/(x)] X Var[^5/(x)] = i. (17) 

The variance of the Gaussian beam along the propagation 
for any angle a is given by the equation [T3l 

M cr^ cos^ a sin^ a 
Var[J-«/(x)] = + (18) 

Let's consider the relation which gives the generalized 
incertitude principle from two applications of the opera- 
tor on f(x) (with first J^a and then J^^) in the case of a 
Gaussian function and displayed in Fig [Sb : 

r./M r^./M cT^cos^tt sin^ a 
Var[/(x)] X Var[J',/(x)] = + (19) 

and compare to that given by Shen [isj in the general 
case of any function (f): 

Var[(/)] X Vd^TlTa^] > (20) 



The generalized uncertainty principle for a rectangular 
function 

The rectangular function s{x) is a peculiar case. Be- 
cause of the abrupt discontinuity, the variance of J^a^i^) 
is infinite whatever a ^ 0. The intensity profile in Fig[T] 
never vanishes completely, from the Fresnel to the Fraun- 
hofer regime [l^. In terms of quantum mechanics, the 
probability for finding the particle anywhere until the 
first moment is not zero. 

To treat diffraction of slit by considering a rectangular 
function is not completely right from an experimental 
point of view, since absorption through the blades of the 
slit induces a not abrupt truncation of the incident wave 
front [sl and thus a finite variance. To take into account 
this effect, we could apply the FrFT to a rectangular 
function convoluted by a Gaussian function g{x)^ 

Ta[s{x) g{x)], (21) 

which would lead to intensity profiles with finite exten- 
sion. Too large Gaussian functions smooth the diffraction 
pattern and reduce the wave extension but may make 
the dark spot disappear. It is difficult to sufficiently re- 
duce the wave extension without vanishing the dark spot, 
which is observed experimentally. We thus measure the 
variance over a limited area centered at the maximum 
intensity. For each a, the full width at half maximum 
is measured and the rms variance is calculated over 2.35 
times the FWHM. This is justify from an experimental 
point of view since, in most cases, the lack of intensity 
or finite sizes of detectors do not allow us to measure 
the diffraction pattern of slit far from the direct beam. 
The result is displayed in Fig [5]i) . A clear minimum is 
obtained for a = /3 = which shows that a successive 
coherent diffraction of two apertures focuses beams in 
the Fresnel regime. This is the basic idea of Fresnel zone 
plates used in X-rays. 



Fractional Fourier transform and diffraction of a 
periodic modulation 

Let's consider a 2D periodic modulation defined by a 
single wave vector go = such as: 

p{x,y) = pocos{qox), (22) 

which gives rise to two Bragg reflections at ig'o by diffrac- 
tion. By continuously varying a from to the FrFT 
simply allows us to calculate the continuous evolution of 
the diffraction pattern from the real space to the recip- 
rocal space (see Fig[6j). 
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FIG. 6: a) 2D periodic modulation defined by a ID wave vec- 
tor qo — '^x. b) Corresponding Fractional Fourier transform 
displayed for a g]0, In the reciprocal space {a = the 
two Bragg reflections are located at g = db^o- The fringes 
around q = ±qo are due to the finite size object. 

Successive diffraction of two objects 

As discussed in the introduction, the use of rectan- 
gular slits is necessary to obtain coherent X-ray beams 
from synchrotron sources and their location relative to 
the diffracted object may influence diffraction patterns. 
To quantify this effect, the diffraction of two successive 
objects has to be taken into account: diffraction of a rect- 
angular function s{x) followed by the modulation 
as respect to (a, /3) G [0, 

J'aiJ'M^)] XP(^)]- (23) 

If the sample is located in the Fraunhofer regime of the 
aperture {a = 7r/2) and the detector in the Fraunhofer 
regime of the sample (/3 = 7r/2): 

[T^ [s{x)] X p{x)] = s{x)® [p{x)] . (24) 

The Bragg peak will mainly display the diffraction pat- 
tern of the rectangular function in the Fresnel regime. 
If now the sample is located in the Fresnel regime of 
the aperture (/3 ~ 0) and the detector in the Fraunhofer 
regime {a = ^), eg. [231 gives: 

J-^ [To[s{x)] X p{x)] = {s(x)\ ® [p{x)]. (25) 

In that case, the reflection profile is a convolution of the 
FT of the aperture with the FT of the periodic modula- 
tion. The profile versus the distance between the aper- 
ture and the sample is summarized Fig [71 It is worth- 
while to note that the Fig. 7d corresponds to the inverse 
of the Fig [TJ By observing Bragg reflection, the double 
diffraction is similar to a time reversal operator. 

The product of variances Yd.T[Fas{x)] times 
\d.TF2L[J^oc[s{x)] X p{x)] is displayed in Fig [H in the 
case of a Gaussian beam and a rectangular aperture. It 
is clear that to increase the width of Bragg reflection, 
the slit has to be as close as possible from the sample. 

To conclude, the fractional Fourier transform ap- 
pears to be appropriate to treat coherent diffraction. 




FIG. 7: a) and b) Successive diffraction of a rectangular func- 
tion and a periodic modulation versus the slit-sample dis- 
tance /3) and the sample-detector distance a), c) for 
a = 0, the diffraction pattern corresponds to the Fig. 5c 
times cos(goa;) d) For a = 7r/2, the Bragg reflection at q=go 
displays a cardinal sinus squared profile for /3 = and a Fres- 
nel diffraction profile when /3 = 7i/2. 




FIG. 8: Product of variances in the case of the diffractions 
two objects: YaT{Ta[T^[s{x)] x p{x)]) with {a, /3) G [0, ^] in 
the case of a) a Gaussian incident beam or b) in the case of 
a rectangular function. The rms variance is calculated from 
the Bragg peak maximum. With = 2q^a^ (1 + exp[^^(7^]) 
et 6 = 2/(^V') (l + expfeV^]). 



Especially, the property of continuity of the FrFT 
could be useful for iterative reconstruction algorithms as 
ptychogr aphy [l5| . 
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